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ABSTRACT 

We present a detailed study of the Faraday depth structure of four bright (> 1 Jy), 
strongly polarized, unresolved, radio-loud quasars. The Australia Telescope Compact 
Array (ATCA) was used to observe these sources with 2 GHz of instantaneous band- 
width from 1.1 to 3.1 GHz. This allowed us to spectrally resolve the polarization 
structure of spatially unresolved radio sources, and by fitting various Faraday rota- 
tion models to the data, we conclusively demonstrate that two of the sources cannot be 
described by a simple rotation measure (RM) component modified by depolarization 
from a foreground Faraday screen. Our results have important implications for using 
background extragalactic radio sources as probes of the Galactic and intergalactic 
magneto-ionic media as we show how RM estimations from narrow-bandwidth ob- 
servations can give erroneous results in the presence of multiple interfering Faraday 
components. We postulate that the additional RM components arise from polarized 
structure in the compact inner regions of the radio source itself and not from polarized 
emission from Galactic or intergalactic foreground regions. We further suggest that 
this may contribute significantly to any RM time-variability seen in RM studies on 
these angular scales. Follow-up, high-sensitivity VLBI observations of these sources 
will directly test our predictions. 
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1 INTRODUCTION 

Radio-loud Active Galactic Nuclei (AGN) eject powerful 
jets of relativistic plasma whose polarized, non-thermal 
synchrotron radiation can be used as a probe of the 
magneto-ionic material along the entire line of sight be- 
tween us and the source of emission. Many studies have 
used these extragalactic background sources to study the 
stre ngth and structure of magnetic fields in our Galax;y 



(e.g. 'Brown ct al.' 2007VTaylor ct al!''2009"; Mao ot al.l l201C 



■ Liaiaxy 
al] |2010l : 



Van E ck ot al . 2011; Ha rycY^LSniith^gt al. 2011 ), other galajc- 
ies (e .g. ICaensler et all booi T iMao et al.l liooa iFeain et al.l 
20091) and in galax y clusters (e.g. iLaing et al.l 1200^ 
Bonafede et~al] l2010l : iPizzo et al.l boill ). Future studies 
on new revolutionary instruments such as the Australian 



Square Kilometre Array Pathfinder (ASKAP) and the 
Square Kilometre Array (SKA) will rely on these back- 
ground sources to probe the strength, structure and evo- 
lution of cosmic magnetism in unpre cedented detail (e.g. 
iBeck fc Gaensleill2004l : lGaenslerll2009l ). 



In this paper, we present a detailed study of the polar- 
ization and rotation measure (RM) properties of four bright, 
unresolved, strongly polarized, radio-loud AGN. Using the 
new Compact Arra y Broadband Backend (CABB) system 
(|Wilson et all 1201 ll ) on the Australia Telescope Compact 
Array (ATCA), spectropolarimetric studies of these AGN 
were performed using 2 GHz of instantaneous bandwidth on 
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the upgraded receiver system from 1.1 to 3.1 Gh£]. All-sky 
RM surveys such as the planned Polarization Sky Survey of 
the Universe's Magnetism (POSSUM) on ASKAP will mea- 
sure the RM s of ^ 3 million extra galactic radio sources over 
30,000 deg^ (|Gaensler et al.ll2010h . POSSUM will likely have 
300 MHz of instantaneous bandwidth covering the frequency 
range from 1130 to 1430 MHz. Proper interpretation of the 
results from this huge dataset will require extensive testing 
of the algorithms used to accurately extract the polariza- 
tion and RM properties of individual sources. The ATCA 
is an ideal instrument for this process, whilst also providing 
new and unique insights into the Faraday depth structure 
of extragalactic sources due to its wide-bandwidth and high 

spectral resol ution. 

Following ISokoloff et al.l (| 19981 ). and references therein, 
we define the complex linear polarization as 

P = Q + iU =ple^'^ (1) 

where I, Q, U are the measured Stokes parameters and 'I' 
is the observed polariza tion angle. We use the notation of 
iFarnsworth et al I (|201ll ) in defining q = Q/I and u ^ U/I, 
so that the measured magnitude of the degree of linear po- 
larization is 

p = + ^2 (2) 
and the polarization angle is 

1 u 

5* = - arctan — (3) 

2 q 

Taking the fractional values decouples depolarization effects 
from simple spectral index effects in analysing the depen- 
dence of polarization with wavelength. It also minimises er- 
rors in the estimate of the RM usin g the RM synthesis tech- 
nique (|Brentiens fc de Bruvn|[2005l ). 

The observed polarization angle is modified from its 
intrinsic value (^I'o) by the effect of Faraday rotation, caused 
by magneto-ionic material between the source of polarized 
emission and the telescope. If there are different regions of 
polarized emission sampled within a single resolution ele- 
ment then each of these regions will likely experience dif- 
ferent amounts of Faraday rotation. Hence, to describe the 
Faraday rotation of a particular region of polarized emission 
we use the Faraday depth (Burn 1966) 

/■telescope 

(j) = 0.81 / nB-dl rad m"^ (4) 

^ emission 

where n is the free electron density (in units of cm~^), B is 
the magnetic field (in /jG) and I is the distance alon g the line 
of sight (in parsecs). iBrentiens fc de Bruvnl (|2005l ) define a 
"Faraday thin" source as one in which A^A<^ ^ 1, and a 
"Faraday thick" source in cases where A^A<jf> ^ 1 (where 
/S.(f) is the extent of the source in Faraday depth and A is the 
wavelength) . 

In the simplest possible scenario, in which there is a 
background source of emission and only pure rotation due 
to a foreground magneto-ionic medium, then the Faraday 
depth is equal to the RM and we get 

* = *o + RMA^ (5) 
^ http:/ /www. atnf.csiro.au/observers/memos/AT39. 3_128.pdf 



Depolarization from radio sources, where the degree of 
polarization decreases with increasing wavelength, is typ- 
ically modelled as a si ngle RM compone nt with externa l 
Faraday dispersion fe.g. iTribble 1991; Rossetti et al.ll2008l '). 
Our new wide-bandwidth data allow a detailed investigation 
of the case of multiple interfering RM components either 
along the line of sight or intrinsic to the source itself. Mul- 
tiple RM components can cause both increases and/or de- 
creases inp(A'^) with A^ as well as, bu t not a lways, deviations 
from a linear 'I'(A^) behaviour. ISlvshI (|l965l ) first employed a 
two component mo del to explain pola rization measurements 
of Cygnus A while [Goldstein fc Reedl (l984) applied a simi- 
lar model to 3C 27. A more recent study bv .Law et al.l (|201ll ) 
showed that multiple RM components could be identified in 
extragalactic point sources using the RM s ynthesis technique 
on wide-band data from 1.0 to 2.0 GHz. IFarnsworth et ahl 
(2011) highlighted the importance of describing both ^{}?) 
and p(A^) in determining the correct Faraday depth struc- 
ture of extragalactic sources using a combination of data at 
350 MHz and 1.4 GHz. Following on from this work, we con- 
clusively show the effect of multiple RM components in two 
extragalactic sources by considering several different Fara- 
day rotation models to simultaneously describe both \&(A^) 
and p(A^). 

In Section 2, we describe the observations, source se- 
lection and calibration process. Section 3 outlines the RM 
synthesis technique while Section 4 describes the various po- 
larization models employed as well as our method for dis- 
criminating between models. Section 5 presents our results 
for each source in order of increasing RM complexity. Sec- 
tion 6 discusses the implications of this work and we list 
our conclusions in Section 7. Throughout this paper, we as- 
sume a cosmology with Ho — 71 km s~^ Mpc~^, Q,m ~ 0.27 
and nA = 0.73, and define the spectral index, a, such that 
the observed flux density (7) at frequency y follows the re- 
lation Ii, oc i/^". 



2 OBSERVATIONS AND DATA REDUCTION 

The sources presented in this paper (Table [1]) were observed 
on 9 Jan 2011 and 20 Jan 2011 with the ATCA from 1.1- 
3.1 GHz (with 1 MHz spectral resolution) in the 6A array 
configuration as part of a larger project to find suitable po- 
larization calibrator sources for the ASKAP telescope. Fig- 
ure[T]shows the typical iiii-coverage for sources in our exper- 
iment. The ATCA has six 22m antennas with linear feeds 
and a maximum baseline of 6 km providing an angular reso- 
lution of ~10" at 1.4 GHz. It has recently undergone a major 
upgrade with the c apabilities of the new CABB system de- 
scribed in detail bv I Wilson et al.l (|201ll '). 

Since there are no 1.4 GHz polarization surveys 
of the southern s ky at a resolution better than 36' 
(jTestori et al.|[200^ . we compiled a list of sources suitable 
for ASKAP polarization calibration from archival obser- 
vations. Candidate sources were selected by searching the 
lATCA online archiveQ for sources that had Stokes I flux 
densit ies greater than 1 J y at 843 MHz in the SUMSS cata- 
logue (|Mauch et al.ll2003l ). The extracted sources, which had 

^ http://atoa.atnf.csiro.au/ 
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Table 1. Observed sources (listed in order of increasing RM complexity) 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


Source 


RA 


DEC 


I 


b 


beam 


pa 


z 


Date 


^obs 


Ao 




a 


Name 


[J2000] 


[J2000] 


n 


n 


r X "] 


[°] 






[min] 


[m] 


[mjy] 




PKS B1903-802 


19:12:40.0 


-80:10:05.9 


314.0 


-27.6 


11 X 9 


-10 


0.50 


2011-01-20 


23 


0.141 


1104 


-0.04 


PKS B0454-810 


04:50:05.4 


-81:01:02.2 


293.9 


-31.4 


12 X 9 


+60 


0.44 


2011-01-20 


16 


0.108 


1009 


+0.46 


PKS B1610-771 


16:17:49.2 


-77:17:18.5 


313.4 


-18.9 


10 X 9 


+30 


1.71 


2011-01-20 


12 


0.137 


3022 


-0.36 


PKS B1039-47 


10:41:44.6 


-47:40:00.1 


281.4 


+9.7 


12 X 9 





2.59 


2011-01-09 


22 


0.164 


1699 


-0.38 


Column designation: 


1 - Source name 


(lAU B1950.0); 2 - Right Ascension 


in J2000 


coordinates; 3 - 


Declination in J2000 coordinates; 4 


Galactic 


longitude 



in degrees; 5 - Galactic latitude in degrees; 6 - Synthesised beam in arcseconds; 7 - Synthesised beam position angle in degrees; 8 - redshift, taken from 
IStickel et al and references tlierein, except for PKS B1039-47 {O. Titov, private communication); 9 - Date of observations; 10 - Total integration 

time in minutes; 11 - Weiglited mean A; 12 - Total intensity at Aq ; 13 - Spectral index (a), defined as oc A~", calculated from our observations. 



been observed using the old narrow-band 128 MHz system, 
were then c alibrated and im aged in Miriad using standard 
techniques (ISault et al.ll 19950 . Eight of the brightest sources 
found in polarized intensity were selected along with two un- 
polarized sources (fractional polarization < 0.1%) for high 
precision polarization observations with the new wide-band 
CABB system on the ATCA. We selected four sources for 
detailed analysis in this paper that had reliable calibration 
across the full 2 GHz band and displayed a range of RM 
complexity from single to multiple RM components. 

In order to avoid any frequency-dependent calibration 
efTects, the data for each source were split up at 128 MHz 
intervals and calibrated using the Miriad software pack- 
age. For the absolute flux density scale correction and band- 
pass calibration, we used a single observation of the ATCA 
primary flux calibrator PKS B1934-638 on each day. Flag- 
ging was done in an autom ated fashion (after bandpass cal- 
ibration) using MiRFLAG ('Middclberd bOOd '). Some minor 
manual flagging was sometimes required afterwards. In total 
~10— 30% of the 2 GHz band was lost due to radio- frequency 
interference (RFI), mainly between 1.1 and 1.7 GHz. The 
leakage and complex gain solutions were determined for each 
individual 128 MHz sub-band before recombining the entire 
band. 

One of our targets, PKS B1903-802, was used on 20 
Jan 2011 as the secondary calibrator to correct for any at- 
mospheric phase variations as well as the polarization leak- 
ages (21 X 1 min cuts, with one cut approximately once ev- 
ery hour, covering 340° of parallactic angle). Because PKS 
B1903-802 is strongly polarized across the entire 2 GHz 
band, the XY phase variations on the reference antenna 
could be determined and we were able to calibrate the ab- 
solute polarization position angle to within ±1°. These so- 
lutions were then copied to the other sources and amplitude 
and phase self-calibration was performed on PKS B0454-801 
and PKS B1610-771. Our other target, PKS B1039-47, was 
used for polarization calibration on 9 Jan 2011 (11 x 2 min 
cuts every hour, with a parallactic angle coverage of 210°). 

Our calibration strategy allowed us to calculate the on- 
axis leakages in 128 MHz intervals. The results show that 
the real part of the leakage solution, which gives informa- 
tion about how the linear feeds deviate from perfect orthog- 
onality (feed misalignment), is constant to better than 0.5% 
across the entire band. The imaginary part, which probes 
the feed ellipticity (i.e. how the feeds differ from perfectly 
linear feeds), has maximum deviations of up to ~ 2% be- 
tween 1.2 and 1.8 GHz. We find little variation between the 



leakage solutions on different days suggesting that the polar- 
ization performance of the wide-bandwidth system is stable 
on timescales of at least a few weeks. 



3 EXTRACTING THE POLARIZED SIGNAL 

We first created uniformly-weighted /, Q and U images for 
each source in 10 MHz intervals and then de convolved thes e 
maps using the Hogbom CLEAN algorithm l|Hogbomlll974l ') . 
In order to avoid any resolution dependent efTects across the 
2 GHz band, we smoothed each image to the resolution at 
the lowest frequency (see Table [T]). Since all sources were 
spatially unresolved, we took the emission at the position of 
the source in the Stokes I image and created a table of q 
and u as a function of A^, at intervals in A'^ corresponding 
to 10 MHz steps. Errors in each channel measurement were 
assigned using the rms noise from a small area around the 
source position in the clean-residual images. 

We then used the RM synthesis technique 
(iBrentiens fc de Bruvnl |2005| ) to extract the polarized 
signal over a wide range of possible Faraday depths. 
RM synthesis is a powerful analysis tool for polarimetry 
data since it helps overcome problems such as bandwidth 
depolarization and nn polarization angle ambiguities. A 
spectrum of complex polarization versus Faraday depth was 
created from these data using the equation 

JV JV 

F(0) = ^^,P,e-^^^(^-^oy^^, (6) 
j=i j=i 

where N is the number of input maps, Pj is the complex 
polarization at channel j and Wj are the weights (inverse 
square of the rms noise). Our reference wavelength (Ao) is 
defined as 

JV JV 

A^ = ^»,A?/^»,. (7) 

j=i j=i 

Essentially this derotates the q and u data for a particular 
assumed RM value and then sums the signal across the band; 
at the correct RM, the channels add coherently giving the 
maximum polarized intensity and the sensitivity of the full 
bandwidth. 

Figure [2] shows the Rotation Measure Spread Function 
(RMSF) for the observations on both days. The RMSF is 
the normalised response function in Faraday depth space 
to the incomplete A'^-sampling (i.e. with perfect A'^ cover- 
age this would be a delta function). More data were flagged 
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Table 2. RM Synthesis Capabilities of 1.1-3.1 GHz ATCA band 



(1) (2) (3) (4) (5) 

Frequency Resolution 50 max-scale |0maa;| 

[GHz] [MHz] [rad m'^] [rad m-^] [rad m-2] 

1.1-3.1 1.0 60 340 13000 



undergoes a different amount of Faraday rotation compared 
to the polarized emission coming from the near side, causing 
depolarization when summed over the entire region. For the 
particular case of a uniform slab we have 

where R is the Faraday depth through the region. 

2. Internal Faraday dispersion (IFD): this occurs when 
the emitting and rotating regions also contain a turbulent 
magnetic field. In this case depolarization occurs because the 
plane of polarization experiences a random walk through the 
region. For identical distributions of all the constituents of 
the magneto-ionic medium along the line of sight, it can be 
described by 

where, in this case, ^'o = n/2 for a purely random 
anisotropic magnetic field and <;rm is the internal Faraday 
dispersion of the random field. 

3. External Faraday dispersion/beam depolarization: 
this occurs in a purely external, non-emitting Faraday 
screen. In the case of turbulent magnetic fields, depolariza- 
tion occurs when many turbulent cells are within the syn- 
thesised telescope beartjfj- On the other hand, for a regular 
magnetic field, any variation in the strength or direction of 
the field within the observing beam will lead to depolariza- 
tion. Both effects can be described by 

where ctrm is the dispersion about the mean RM across the 
source on the sky. 

A fourth possibility is a changing degree of polarization 
due to multiple interfering RM components, either along the 
line of sight or on the plane of the sky on scales smaller than 
our spatial resolution. 

To model our data, we simultaneously fit both the 
g(A^) and u(A^) data to the different polarization models 
listed above. We first tried a simple one-component RM 
model, which cannot describe any variation in the degree 
of polarization. We then tried to account for changes in 
p(A^) by fitting a single RM component plus depolariza- 
tion model, and in cases where this could not adequately 
describe the data, we then tried multiple RM-component 
models. In order to limit the number of models to investi- 
gate, we mainly considered models of either solely Faraday 
thin components (e.g. emission from the radio galaxy only) 
or models of one Faraday thick component (e.g. Galactic 
slab or mixed emitting and rotating region in the source) 
plus Faraday thin component (s). We believe that these mod- 
els are the most physically reasonable cases for spatially 
unresolved extragalactic radio sources. However, we do not 
consider differing spectral indices of individual components. 
Hence, multiple component models are simply constructed 
asP = Pi+P2 + --- + PN. 



However, see iTribbld lll99ll ) for a detailed description of what 
happens when this assumption does not hold. 



Column designations: 1 - Instantaneous frequency coverage in GHz; 2 - 
Spectral resolution in MHz; 3 - Resolution in Faraday depth space; 4 - 
maximum detectable Faraday thickness; 5 - Maximum detectable Faraday 
depth. 

on 9 Jan 2011 than on 20 Jan 2011, so that the RMSF 
for 9 Jan has slightly stronger sidelobes than on 20 Jan. 
Table [2] lists, for our observations, the Faraday depth res- 
olution (50 ^ 2\/3/A(A2), i.e. the FWHM of the RMSF), 
the largest detectable scale in Faraday depth space (max- 
scale ~ ■""-^mL) '^^'^ maximum observable Faraday depth 
(^max ~ \/3/<5A^), where A(A^) is the total bandwidth in 
A^-space (i.e. A(A^) = A^ax ~ -^min) and is the channel 
width. 

For each source, we initially searched for polarized 
power from -ii4>max at Faraday depth intervals of 10 rad m~^ 
which corresponds to ~6 Faraday depth intervals per 54> (see 
Table |2]). No significant power was found at large values 
of 101. For the rest of our analysis we restricted our range 
to ±1000 rad m~^ at 1 rad m~^ intervals (~ 60 Faraday 
de pth intervals per 5 (j)). We use RMCLEAN, as described 
bv iHeald (|2009l '). to deconvolve the "dirty" RM spec- 

trum in an attempt to recover information lost due to the 
incomplete frequency coverage. For each source, we cleaned 
down to the rms noise level {aq,u) listed in Table 3. In prin- 
ciple, the RM synthesis technique should be able to detect 
regions that are extended in Faraday depth space as long as 
the region extends beyond the FWHM of the RMSF. 



4 MODELLING PROCEDURE 

To model the polarized signal in the presence of Faraday 
rotation in the simplest case, we use the equation 

where po is the intrinsic degree of polarization of the syn- 
chrotron emission, \l/o is the intrinsic polarization angle at 
the source of the emission and the RM describes the Faraday 
rotation caused by the foreground magneto-ionic material 
with the sign indicating whether the line-of-sight magnetic 
field is pointing towards us (positive RM) or away from us 
(negative RM). The data for all our sources show changes 
in the degree of polarization across the observed wavelength 
range so we now consider the possible mechanisms behind 
this effect. 

Depolarization towards longer wavelengths can occur 
due to mixing of the emitting and rotating media, as well as 
from the finite spatial resolution of our observations. There 
are three commonly listed depolarization mechanisms (see 
ISokoloff et ai] (Il998l ) for more detailed discussion of each 
case) . 

1. Differential Faraday rotation (DFR): this occurs 
when the emitting and rotating regions are co-spatial and 
are in the presence of a regular magnetic field. The polar- 
ization plane of the emission at the far side of the region 
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4.1 Model-fit Evaluation 

We utilised the maximum likelihood method to find the best- 
fit model parameters. The results from RM synthesis were 
used as a guide to our initial guesses for the RM and frac- 
tional polarization values. Each data point in the fit was 
weighted by the inverse square of the rms noise from the 
clean-residual image. 

The likelihood is the probability of obtaining the data, 
d, given a model of the source and some characterisation of 
the noise. Our data in this context are g(A^) and w(A^). For 
example, in the one-component model, we adopt qmodei,i ~ 
po cos(2«'o + 2RMAf ) and u„,odei.i = Po sin(2*o + 2RMA?), 
and we assume that qi = qmodeis + rii, where Ui is Gaus- 
sian noise for channel i. The prior likelihood of a particular 
RM value for an observation of a single channel i under the 
assumption of Gaussian noise is 



1 f {q^ qmodel.i) (ui ^modeZ.i) 

exp ' 



(12) 

where Oq^u is the single channel rms. If we have a total of A'' 
channels, the prior likelihood is now 



P(d!RM) ^\\Pi{d,\RM) 



(13) 



We used the Mathematic^ function Nonlinearmod- 
ELFIT to find the maximum of Egn.fTSl L = max{P{d\RM)). 
In the case of multiple models giving good fits to the data, 
we then used the Bayesian Information Criterion (BIG) to 
distinguish the goodness-of-fit between different models with 
different degrees-of- freedom ()Schwarz.,1978, : ,Trottai2008i ') : 

BIG = -21ogP(d|fe) « -21ogL + felogAr (14) 

where k is the number of free parameters in the model. 
Hence, models with more parameters are heavily pe- 
nalised given the large number of data points. We con- 
sider BIGmodcii ~ BIGmodci2 > to siguificautly favour 
model 2 (at 99% level) across the measured parameter space. 
Bayesian model comparisons require an alternative model 
against which the comparison is made (i.e. a model cannot 
be rejected unless an alternative explanation is available that 
better fits the observations). In order to give a quantitative 
measure of how each individual model fits the data, we also 
calculate the reduced chi-square (xl) goodness-of-fit values, 
where xt is obtained by dividing the sum of squared resid- 
uals by (Tg „ and the number of degrees of freedom. 



5 RESULTS 

We present the sources in order of RM complexity (as listed 
in Table [T]) with PKS B1903-802 having the simplest RM 
structure and PKS B1039-47 having the most complex. All 
models with their best-fit parameters for each source are 
listed in Table 3 (with the most-likely model highlighted in 



* Wolfram Research, Inc., Mathematica, Version 7.0, Champaign, 
IL (2008). 



bold). The associated errors of each parameter are formal 
fitting errors which are calculated from the square root of 
the estimated error variance of each parameter. These errors 
have little meaning when the incorrect model is applied. 
Note that all polarized components listed are found with 
high significance. For example, the weakest polarized model 
component of 0.6% listed in Table 3 has a signal to noise 
ratio (SNR) of ~100. 

In Figure [31 we show the results of RM synthesis and 
RMGLEAN for all four sources. For each source we list 
the RM at the peak degree of polarization with its associ- 
ated error. The uncertainty in the peak RM is calculated 
as the FWHM of the RMS F divided by twice the SNR 
(|Brentiens fc de Bruvnl l2005i ) . So for example, if we have 
an SNR of 600 and an RM resolution of 66 rad m~^, then 
the quoted uncertainty i n the peak RM is 0.06 rad m~^. 
However, as described bv iLaw et al.1 (|201ll ). the accuracy of 
any individual RM-component value cannot be specified to 
better than the RM resolution due to the uncertainty in the 
distribution of components within the RM beam. The mean- 
weighted reference wavelength (Ao) changes from source to 
source mainly due to different amounts of flagged data across 
the band for each source. 

Faraday depth spectra for both PKS B1903-802 
(Fig.[3ji) and PKS B0454-810 (Fig.|3}3) appear to be broadly 
consistent with a single RM component. The asymmetric 
distribution about the peak for PKS B1610-771, as well as 
the distribution of clean-components, indicates the presence 
of more than one RM component (Fig.[3j:). PKS B1039-47 
has a distinct secondary peak at ~100 rad m^^ and the 
clean-component locations suggest the presence of three or 
more RM components (Fig. |3ji) . We now discuss the Fara- 
day rotation model-fits to the (?(A^), u{\'^) data for each 
source which are completely independent of the RMGLEAN 
results. 



5.1 PKS B1903-802 

A simple RM fit (Eqn. [S]), as shown in Figure UJ provides 
a reasonable description of the ^'(A'^) data but the p vs. A'^ 
data clearly deviates from a constant degree of polarization. 
An external Faraday dispersion model (Eqn. Illjl provides 
an excellent fit to data, with our best fit model, shown in 
Figure [5] giving a polarized intensity of 5.14 ± 0.04%, with 
a foreground RM of -fl8.1 ± 0.1 rad m~^ and an external 
dispersion in RM across the source of4.7±0.1 rad m'^ How- 
ever, in Figure|6]we can see that a two RM-component model 
also provides a very good description of the data with an RM 
of -f 16.1±0.3 rad m~^ for the stronger polarized component 
(~4.5%) and an RM of 39.3±2.0 rad m"^ for the second po- 
la rized component ( ^0.6% ). This supports the conclusion 
of iFarnsworth et al.l l|201ll ') that modelling of both polariza- 
tion amplitude and polarization angle is required in studies 
of Faraday rotation. 

In this case, neither the BIG nor the xt help us to clearly 
discriminate between models. The two models do not dif- 
fer significantly over the measured parameter space so we 
cannot state with confidence which one is correct, although 
we favour the simpler external Faraday dispersion model. 
Lower frequency observations, from 700 MHz to 1 GHz on 
ASKAP for example, would clearly discriminate between the 
two models since for the two RM-component model there is 
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a departure from a linear 'I' vs. A'^ relationship over this 
range. The difference can be seen most clearly in the g(A^) 
vs. u(A^) plots (i.e. compare Fig. [5] and Fig. [6]). 

5.2 PKS B0454-810 

As can be seen in Figures [7] & [HI both single-component RM 
models (with and without a depolarizing screen) provide 
poor fits to the data. They both determine approximately 
the same RM (+37.8 ± 0.2 rad m~^) but cannot explain 
the observed decrease in the degree of polarization towards 
the shortest wavelengths. The two-component model, shown 
in Figure O does much better at describing the p(A^) data 
while also providing a good description of ^^(A^). We note 
that the RM of the strongest polarized component is now 
significantly different (+29.2±1.1 rad m"'^) from that found 
in the one-component models or from the peak in Faraday 
depth inferred from RM synthesis in Figure |3Ib). 

While both the reduced chi-squared and BIG values 
strongly favour the two RM-component model over the sin- 
gle component models (Table 3), on closer inspection it is 
clear that the p(A^) data at the shortest wavelengths ob- 
served are not very well fit by the two RM-component model 
either (Fig.|5}. Hence, we conclude that the data are not well 
fit by any of the three models listed in Table 3. 

We consider a possible alternative explanation for this 
source in terms of polarization propagation effects as a func- 
tion of optical depth within the source. Specifically, the case 
of one optically thick and one optically thin component; the 
optically thin and strongly polarized component would dom- 
inate the emission at the longer wavelengths (A^ > 0.025 m^) 
while the observed depolarization could be explained by 
external Faraday dispersion. At the shorter wavelengths 
(A^ < 0.025 m'^), the optically thick and weakly polarized 
component would become more dominant coupled with the 
optically thin component becoming fainter leading to a de- 
crease of the observed degree of polarization in a non-trivial 
manner. This idea is supported by the inverted spectrum 
indicative of a synchrotron self-absorbed region caused by 
either multiple, discrete spectral components or a smooth 
distribution of magnetic field and electron density along the 
jet. If this model were correct, then we would expect at even 
shorter wavelengths, as the emission spectrum turns over, 
the degree of polarization should increase again. However, 
detailed modelling of the polarized radiative transfer from 
such a region is required for a quantitative analysis and we 
defer such a study for a later paper. 

5.3 PKS B1610-771 

Figure [TU] shows how a simple one-component model pro- 
vides a very poor fit to the complex polarization data for 
this source. A depolarizing screen model, shown in Figure fTTI 
gives a good fit at the short wavelengths but fails to ade- 
quately fit the long wavelength data. These plots further 
highlight two interesting features of the data. First, it is 
clear that the slope of the >I'(A^) relationship gets steeper 
towards longer wavelengths and, second, while thep(A'^) plot 
shows that emission is strongly depolarized it also shows ev- 
idence for a reversed trend of increasing polarisation with 
A^ at longer wavelengths. Any simple depolarization model 
cannot explain both these effects. 



A two RM-component model, shown in Figure [121 pro- 
vides a much better description of the data as demonstrated 
by the corresponding value of 1.04. The fit accounts 
for both the changing slope of ^(A^) as well as the in- 
creasing p(A^) for A^ > 0.05 m^ The best-fit RMs are 
-1-107.1 ± 0.2 rad m"^ and -1-78.7 ± 0.4 rad m"^ for the 
first and second component, respectively. The BIG strongly 
favours the two RM-component model over the depolarizing 
screen model (see Table 3 for values). 



5.4 PKS B1039-47 

This is the most striking source in terms of complex po- 
larization structure. Figure [131 shows how both ^'(A^) and 
p{y^) display non-linear, oscillatory behaviour indicative of 
multiple RM components. The RMGLEAN spectrum also 
indicates the presence of multiple RM components (Fig.[3}l). 

We list a single RM-component model in Table 3 for 
completeness but do not show the fit. A two RM-component 
model (both components Faraday thin) also provides a poor 
description of the data (Fig. I13[) with a reduced-x^ value 
of 14, and this can be seen quite clearly in the p(A^) dis- 
tribution. We then tried models with one Faraday thin and 
one Faraday thick component where the Faraday thick com- 
ponent was described by either Eqn.[^or Fan. 1101 In both 
cases, neither model converged to an acceptable solution and 
they are not shown here. 

Our next approach was to fit models with three RM 
components, which in the case of all Faraday thin compo- 
nents have a total of nine parameters (i.e. po, RM, & 
for each component). In order to find a good model, we 
first fixed the parameters of the dominant RM component 
taken from RMGLEAN and let the other six model param- 
eters vary. We then used these results as input guesses for 
the final nine parameter model fit. This returned best-fit 
RMs of -13.1 ± 1.5 rad m'^ -29.8 ± 2.4 rad m"^ and 
-1-68.4 ± 1.6 rad m~^, listed in order of highest to lowest 
polarized fractions for the individual components. This pro- 
vides a good fit, shown in Figure 14, with a value of 1.2 
(Table 3). 

We also tried different combinations of Faraday thin and 
thick components within a three RM-component model with 
the best model shown in Figure 15. Both the BIG and x? val- 
ues are almost identical for both types of three-component 
model listed in Table 3, so they do not help us to clearly 
discriminate between models in this case. More data, at 
longer wavelengths, is required to determine which three- 
component model provides a better fit. A four-component 
model (all Faraday thin) was also tried but did not improve 
the fit. 



6 DISCUSSION 

6.1 Physical origin of the RM components 

We have conclusively shown in the previous section that 
for two sources (PKS B1610-771 & PKS B1039-47) we can 
spectrally resolve multiple polarized components of spa- 
tially unresolved AGN. We now discuss the likelihood of 
these additional RM components coming from regions of 
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(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


Source 


Model 


RMi 


POl 




RM2 


P02 


*02 


RM3 


P03 


*03 






xi 


BIC 






[rad m~^] 


[%] 


[°1 


[rad m — ^] 


[%] 


n 


[rad m~^[ 


l%] 


[°1 


[rad m~"] 


[%] 






PKS B1903-802 


Single RM, no screen 


+ 18.1(1) 


4.98(1) 


-4.4(1) 
















0.006 


1.23 


-3470 




Single RM, foreground screen 


+ 18.1(1) 


5.14(4) 


-4.4(1) 














4.7(1) 




1.03 


-3714 




2 R.M Components, no screen 


+ 16.1(3) 


4.50(8) 


-2.4(3) 


+39(2) 


0.59(8) 


-25(2) 












1.04 


-3779 


PKS B0454-810 


Single R.M, no screen 


+37.9(2) 


3.36(2) 


-48.6(3) 
















0.008 


2.13 


-3250 




Single RM, foreground screen 


+37.8(2) 


3.39(2) 


-48.6(3) 














2.8(6) 




1.60 


-3250 




2 RM Components, no screen 


+29(1) 


2.7(2) 


-33(2) 


+62(3) 


1.1(2) 


+88(5) 












1.07 


-3548 


PKS B1610-771 


Single RM, no screen 


+ 109.0(7) 


3.79(7) 


+69(1) 
















0.005 


97.3 


-2346 




Single RM, foreground screen 


+ 104.3(4) 


5.20(4) 


+73.9(5) 














17.3 (0.3) 




1.41 


-3024 




2 R IVl Components, no screen 


+ 107.1(2) 


3.45(4) 


+83.3(6) 


+78.7(4) 


1.98(4) 


+73(1) 












1.04 


-3845 


PKS B1039-47 


Single RM, no screen 


-12.3(5) 


3.43(5) 


+ 18.8(8) 
















0.005 


127 


-2362 




2 RM Components, no screen 


-9.8(4) 


3.64(4) 


+ 14.3(7) 


+85(2) 


0.68(4) 


0(4) 












14.3 


-2524 




3 H.M Components, no screen 


-13(1) 


3.9(4) 


+31(4) 


-30(2) 


1.7(4) 


0(8) 


+68(2) 


0.7(2) 


+35(3) 






1.23 


-3194 




3 RM Components {1 DFR) 


-13(4) 


2.3(3) 


+37(5) 


+70(1) 


0.7(2) 


+33(3) 


-21(4) 


2.9(1) 


+ 11(15) 






1.25 


-3192 



Column designation: 1 - Source name; 2 - Description of Faraday rotation model used; 3 - RM of first component ; 4 - Degree of polarization of first component; 5 - Intrinsic polarization angle (at A — 0) of 
first component; 6 - RM of second component ; 7 - Degree of polarization of second component; 8 - Intrinsic polarization angle of second component; 9 - RM of third component ; 10 - Degree of polarization of 
third component; 11 - Intrinsic polarization angle of third component; 12 - Dispersion about mean RM; 13 - rms noise level in fractional polarization; 14 - Reduced chi-square goodness-of-fit value; 15 - Bayesian 
information criterion for model comparison. The favoured model for each source is indicated by bold-face font. The error in the final digit of all parameters is indicated by the number in the parentheses. 
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polarized emission along the line of sight or from multi- 
ple polarized regions on the plane of the sky but within 
our synth esised beam. Polarized diffuse Galactic emission 
IXestori et a l. 200§) and polariz ed emission from r adio ha- 
los/relics in galaxy clusters fe.g. iFerrari et al]|2008l ) are the 
most likely candidates for any additional polarized emission 
components along the line of sight. However, both these pos- 
sibilities are highly unlikely for our particular observations 
since first, we do not have the sufficient short itu-spacings 
to detect the smooth Galactic emission and secondly, none 
of the sources studied here have any diffuse X-ray emission 
associated with th em in the ROSAT All Sky Survey (RASS; 
IVoges et al.lll999l '). which effectively rules out the presence 
of any significant emission from galaxy clusters along the 
line of sight. 

If the sources were spatially unresolved, double-lobed 
radio galaxies such that our line of sight to one of the lobes 
travelled through a different magneto-ionic medium, then 
the polarized emission from each lobe would experience dif- 
ferent amounts of Faraday rotation and could show up in 
our data as tw o distinct RM components (e.g. ISlvshin"965l : 
[Goldstein fc Reed f984,'). For blazar type sources we require 
variation in the magneto- ionic medium along the jet because 
we only detect the Doppler boosted emission from the jet ori- 
entated toward us. Many high-resolution studies of such ob- 
jects have shown substantial variatio ns in both RM and po- 
larized intensity on parsec-scales (e.g. Zavala fc Tavloi|2003l . 
I2OO4I : lO'Sullivan fc Gabuzdallioogl : iHovatta et al.ll201ll ). Be- 
low we discuss what is already known for each source and 
why the mostly likely origin of the additional RM compo- 
nents is from the compact inner regions of the radio source. 

PKS B1903-802 is a fiat spectrum radio quasar which, 
from 1.4 to 20 GIfz, maintains a total fiux density of ~1 J y 
while remaining polarized at ~3% (jMurphv et all I2OI0I I. 
From 43 GHz observations with the ATCA, we know that 
this source is unresolved down to at least 0.15" (Rajan 
Chhetri, private communication) which corresponds to a lin- 
ear scale of ~0.9 kpc for the quoted redshift of this source 
(Table [l]). The source is resolved into two total intensity 
components on scales less than 5 mas (30 pc) from observa- 
tions a t 8.4 GHz with t he Australian Long Baseline Array 
(LB A; lOiha et ahllioosl ). No spectral index or polarization 
information is available on these scales but at least one of 
these components must be the origin of the strong polarized 
emission we see in our data. Hence, it is likely that there 
is a significant contribution to the observed RM from the 
immediate environment of the source and its host galaxy as 
well as the Faraday rotation caused by our own Galaxy. 

PKS B0454 -810 is a well s tudied flat spectrum ra- 
dio quasar (e.g. iRicci et all 120041 ') and has an inverted ra- 
dio light curve whic h begins to turn o ver above 100 GHz 
(jBennett et al.ll2003l l. iRicci et al .' (2004') detected polarized 
emission of ~2.6% at 18.5 GHz while MujCBliy ot al. (2010.) 
quote an upper limit of 1.5% from 20 GHz ATCA obser- 
vatio ns. While it has been detected in X-rays ijVoees et al.l 
1 19991). it does not ha ve a 7-ray detection from Fermi 
( Abdo et al1l2009l .[ 201 ih . It is spatially unresolved down to 



at least 5" (~30 kpc) from inspection of the visibilities on 
the longest baselines in our observations. It has also been 
imaged on milliarcsecond scales with both th e VLBI Space 
Observatory Prog ramme ijPodson et al.|[2008l l and the LBA 
showing that it is resolved on scales less 



Ubservatory Jr'rog i 
l|Oiha et al.ll2005l ) 



than 5 mas. iRavner et al] (|200G| ) found the source to be cir- 
cularly polarized, with a > 10a detection at both 1.4 and 
5 GHz, which is indicative of a core-dominated AGN. There- 
fore, all this data supports our assertion that the linearly 
polarized emission we detect is coming from the compact 
inner regions of the AGN jet and provides some weight to 
our explanation for the observed variation in p(A^) being 
due to the combination of a weakly polarized optically-thick 
region and a strongly polarized optically-thin region. 

PKS B1610-771 has been extensively studied across 
a wide range of wavelengths . It is classified a s a fiat- 
spectrum radio quasar (e.g. iHealev et al. 

I l200'if ). IS sig- 

nificantly polarized (1-2% level) at 5, 8 and 20 GHz 
llMassardi et al. 2008[) and is highly polarized (>3%) in the 
optical ( Veron-Cettv fc Veronll2006r). It is co incident with an 
unresolved X-ray source (|Voges et al.lll999l ) and has a GeV 
7-ray detection in both the first and second Fermi-LAT cat- 
alogues ( Abdo et a l. 2009, 2011). No source structure is de- 
tected on scales greater than 0.15" (Chhetri, private comm.), 
which corres ponds to a linear scale of ^1 .3 kpc at its red- 
shift of 1.71 (|Hunstead fc Murdoch|[l980h . PKS B1610-771 
is also seen to exhibit interstellar scintillation at low ra- 
dio frequencies, with a char acteristic time scale of 400 days 
(jCaensler fc Hunsteadll200"ol ). This further supports the con- 
clusion that its fiux is dominated by compact rather than 
extended components. On milliarcsecond scales at 8.4 GHz 
the jet extends in a North- West direction with several bright 
jet knots in total intensity seen out to a projected distance of 
~130 pc (|Oiha et al.ll2010l ). From our analysis of this source, 
we predict that two of these knots are strongly polarized and 
have different RMs. 

PKS B 1039-47 is classified as a flat-spectrum radio 
quasar that is located along a sightline ~10° from the Galac- 
tic plane , and has been meas ured to be ~4% polarized at 
20 GHz (iMassardi et al.ll2008l ). The host galaxy has a mea- 
sured redshift of 2.59 (O. Titov, private communication) and 
there is no X-ray or 7-ray detection for this source. The emis- 
sion structure remains u nresolved down to 0.15" (1.2 kpc), 
and an LBA image fa-om lOiha et~aLl (|2004h shows a jet ex- 
tending to the North- West out to ~20 mas (160 pc), com- 
posed of three bright total intensity regions. Our analysis 
in this paper, which flnds a best-flt three RM-component 
model for this source, suggests that each of these regions is 
polarized. 

Therefore, in the case of PKS B1610-771 and PKS 
B 1039-47, we claim to have detected separate polarized com- 
ponents in the compact inner regions of the jet on parsec- 
scales that are illuminating an inhomogeneous magneto- 
ionic medium in the immediate vicinity of the jet. The 
largest RM difference between different components that 
we have found from our model flts is ~100 rad m~^ (be- 
tween the second and third components in PKS B1039-47). 
This is not inconsistent with recent measurements on mil- 
liarcsecond scales at 1.4 GHz where variations of 10s to 
100s of rad m~^ hav e been measured in several blazars 
(jCoughlan et al.ll201ll '). We can test our predictions for these 
sources through multi-frequency polarization sensitive ob- 
servations with the LBA, which will allow us to spatially 
probe their milliarcsecond scale polarized structure for the 
first time. 
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6.2 RM time-variability 

If we attribute the RM difference between components to 
tlie magneto-ionic material in the immediate vicinity of the 
parsec-scale jet then there are obvious implications for any 
observed time- variability from RM measurements on arcsec- 
ond scales and greater. The polarized and RM structure of 
VLBI jets has b een observed to vary on timescales as shor t 
as months fe.g. iGomez et al.ll201ll : iZavala fc Tavloill200ll ). 
so the relative polarized flux of the individual components 
detected in our ATCA observations may also vary on simi- 
lar timescales. Therefore, the relative difference between the 
individual component RM values may also vary if the polar- 
ized components are moving along the jet and illuminating 
different parts of an inhomogeneous Faraday screen close to 
the AGN. This means that RM time-variability for obser- 
vations on similar angular scales to those presented in this 
paper may be simply due to the observations sampling dif- 
ferent dominant RM components as they move along the jet 
or, if the observations do not have sufficient frequency cov- 
erage and spectral resolution, a complicated combination of 
multiple RM components that does not accurately represent 
any of the com p onent s. 

iLaw et all (|20lj l were able to compare four sources 
from their low spatial resolution, wide-bandwidth 1-2 GHz 
observations with high spatial resolution RM maps from the 
Very Long Baseline Array (VLBA). In general, they did 
not detect the high fractional polarization or high RM val- 
ues seen from 5 to 22 GHz in the VLBA images. This is 
not surprising since the high spatial resolution of VLBA 
observations are less affected by beam depolarization and 
also because RMs on these scales have been observed to 
incre ase with increasing frequency jO'Sullivan fc Gabuzdal 
I2OO9I ). Another possibility is that there may be intermediate 
scale polarized structure that the VLBA is not sensitive to, 
negating the validity of the comparison. Therefore, parsec- 
scale VLBI observations at similar frequency ranges taken as 
close in time as possible to the low spatial-resolution, wide- 
bandwidth observations are required for direct comparison. 
If the results from both these types of observations can be 
linked then it may be possible with multi-epoch polarization 
observations with wide-bandwidth facilities like the ATCA 
to map out the parsec-scale evolution of polarized compo- 
nents as well as the Faraday ro t ating environment in AGN 
jets, as suggested bv lLaw et al.l l|201lf ). 

Even though the majority of the Faraday rotation oc- 
curs as the polarized radiation passes through our Galaxy, 
the RM difference between multiple components as well as 
any observed variability is likely due to the magneto-ionic 
material in the immediate vicinity of the AGN jet. Hence, 
the type of sources studied in this paper may not be suit- 
able as primary polarization calibrators since the values for 
their RMs and degree of polarization are likely to change 
on short timescales. A more stable type of calibrator source 
would be one in which the polarized emission comes from 
extended emission regions such as the lobes which vary on 
much longer timescales. 



6.3 Reliability of RMCLEAN 

For PKS B1903-802 the peak RM found using RMCLEAN 
agrees very well (within ~0.3 rad m~^) with the mean RM 



Table 3. Comparison between narrow-band and wide-band RMs 



Source 




R-Moid ATCA 


RMnvSS 




[rad m~-^] 


[rad m~^] 


[rad m~^] 


PKS B1903-802 


-1-18.1 ±0.1 


+18.2 ± 0.4 


+21.5 ± 1.3 


PKS B0454-810 


+37.8 ±0.2 


+39.9 + 1.2 


+38.9 ± 1.5 


PKS B1610-771 


+107.1 ±0.2 


+134.6 ± 1.0 


+128.5 ± 2.7 


PKS B 1039-47 


-13.1 ± 1.5 


-13.0 ±6.7 


-8.9 ±1.2 


RM„ide: RM of 


main component dc 


:rivGd from 1.1—3.1 


GHz data. 



RMnarrow: RM derived using data from 1304- 1494 MHz. RMnvs 
derived using the same frequency coverage as in lTavlor et all j 20091 



found from our best-fit external Faraday dispersion model. 
In the case of PKS B0454-810, the RM extracted using RM- 
CLEAN differs by ~2 rad m~^ from what we consider to 
be the correct RM for this source. However, the comparison 
in this case is somewhat unreliable since we have not been 
able to conclusively identify the correct polarization model 
for this source. 

There have been some questions raised in the literature 
about the ability of the RMCLEAN method to accurately 
recover multiple RM com ponen ts for sources with c omplex 
Farad ay structure (e.g. Frick et a l. 2010; Farnswort h et al.l 
l201lh . We find that the RMCLEAN method performs poorly 
in recovering the correct RMs for PKS B1610-771 and PKS 
B 1039-47. While it does predict the presence of multiple 
RM components, the clean-component distribution does not 
match what is found through model-fitting g(A^), ii(A^) 
(Fig. [3j;,d). For PKS B1039-47, we investigated whether or 
not the three best-fit model RM components could be recov- 
ered with the same A'^-coverage but with no noise. Figure [TOl 
shows that again, after RMCLEAN, the clean-component 
distribution does not associate polarized power with the cor- 
rect RM model-component locations. On refiection this may 
not be very surprising given that the RM resolution of our 
experiment is ~60 rad m~^ and is therefore unable to re- 
solve multiple RM components which differ by less that this 
value. Thus, in cases of complex RM structure of extragalac- 
tic point sour ces, alternative reconstruction algorithms (e.g. 
iLi et al.ll201ll ) need to be investigated for application in all- 
sky RM surveys. Currently, the best approach for sources 
with complex Faraday depth structure is model-fitting mul- 
tiple RM-component models to the observed ?(A^) and u{X^) 
in order to find the most likely physical model for the source. 



6.4 Some implications for RM surveys 

Taking the same data but restricting the frequency range 
to 20 X 10 MHz channels, centred at 1.4 GHz, allowed us 
to compare the RMs derived from the full 2 GHz band- 
widt h with the previo us narrow-band system on the ATCA 
(e.g.lPeai n et all|200Sl ). Using the fitting procedure employed 
throughout this paper, we find that we get the same RMs 
(within the errors) using 200 MHz of data for both PKS 
B1903-802 and PKS B0454-810. In this case, a single RM- 
component model with external Faraday dispersion still pro- 
vides a better fit over a model without depolarization. In 
Table [3] for each source, we list the RM of the strongest po- 
larized component derived from the full 1.1-3.1 GHz data 
to compare with the RM derived from the restricted range 
of 1.340-1.494 GHz. 
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For PKS B1610-771 we obtain a good fit {xl = 1-2, 
BIC= —475) for a single RM-component with external Fara- 
day dispersion giving a mean RM of +134.6 ± 1.0 rad 
and a dispersion of 10.9±0.5 rad m~^. A marginally poorer 
fit is found using a two RM-component model (x^ = 1-3, 
BIC= -452) with RMi = 4-110.1 ± 1.8 rad m'^ and 
RM2 = -1-79.7 ±2.1 rad m~^. Therefore, in this case we 
would adopt the single RM-component model since the evi- 
dence does not favour the more complex two RM-component 
model. But we already know that the two RM-component 
model is preferred using the full 2 GHz bandwidth data. The 
difference in RM between the single-component model and 
the strongest component in the two RM-component model 
is ~ 25 rad m~^ (Table [3)|. This is quite a dramatic exam- 
ple of how wrong one can be in estimating the RM using 
narrow-bandwidth data. This highlights the important role 
wide-bandwidth data play in determining the correct Fara- 
day depth structure of AGN on these angular scales. The 
best-ft model for PKS B1039-47 using the data from 1.340- 
1.494 GHz has two RM components where the strongest 
component has an RM equal (within the errors) to that 
found from the 1.1-3.1 GHz data. 

Due to our small sample size we are unable to comment 
on how often additional RM components may be detected in 
extraga lactic sources. Fro m a sample of 37 bright, polarized 
sources iLaw et al.l (|201ll ) found that ~25% of sources had 
an extra RM component detected with a significance greater 
than 7cr (~40 mjy). They also found that the polarized flux 
weighted mean RM was simila r to the low resolution RMs 
quoted by iTavlor et al. I (|2009l ). This suggests that sources 
like PKS B1610-771 may be rare where the flux weighted 
mean RM (~97 rad m~^) is significantly different from the 
RM which is derived from narrow-bandwidth observations. 
We also list in Table[3]the RMs calculat ed from our data us - 
ing the same frequency setup as used bv lTavlor et al.l (|2009l ). 
For the simple sources, the RMs agree within 2-3 rad . 

Current and upcoming s pectropolarimetric all -sky RM 
surveys such as GALFACTS (Taylor & Salter"2010') and the 
POSSUM survey on ASKAP (Gaensle r et al. .2Q1Q 1 plan to 
extract RMs from observations with 300 MHz of instanta- 
neous bandwidth near 1.4 GHz. For sources with a single 
RM-component modified by depolarization from external 
Faraday dispersion, observations with 300 MHz of band- 
width will produce the same results as we have found us- 
ing 2 GHz of bandwidth (e.g. PKS B1903-802). In the case 
of sources with multiple RM components, it is strongly rec- 
ommended that modelling of q(A^), u{}^) be undertaken 
instead of using the reconstruction algorithm RMCLEAN. 
Since ASKAP will have much better RFI conditions than at 
the ATCA, we use simulated data from the best-fit models 
for PKS B1610-771 and PKS B1039-47 (highlighted in bold 
in Table 3) instead of the observed data. Using the planned 
POSSUM frequency coverage of 1130-1430 GHz, we find 
that our modelling procedure would recover the correct two 
RM-component model for PKS B1610-771. For PKS B1039- 
47, depending on the quality of the data, it may be diffi- 
cult to determine whether a two or three RM-component 
model provides the best flt. Complementary observations 
from a 300 MHz band at lower frequencies (e.g. 0.7-1.0 GHz) 
would enable us to recover the correct simulated three RM- 
component model in this case. 



I 1039-47 2.1005 GHz 




-60 -40 -20 20 40 60 
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Figure 1. Plot showing the typical MD-coverage for the sources 
presented in this paper (PKS B1039-47 shown). 

7 CONCLUSION 

Using the new wide-bandwidth receivers on the ATCA, we 
have shown that we can spectrally resolve the polarization 
structure of spatially unresolved radio sources. We have 
identified two AGN (PKS B1610-771 and PKS B1039-47) 
where more than one rotation measure (RM) component is 
required to describe the Faraday structure of the source. We 
further demonstrate that modelling of both the polarization 
angle and degree of polarization dependences with wave- 
length squared is essential in determining the true Faraday 
depth structure of extragalactic point sources. We also find 
that the RM synthesis reconstruction algorithm RMCLEAN 
does not recover the correct RMs for sources with multiple 
RM components in our data. 

The most likely origin for the additional RM compo- 
nents in both PKS B1610-771 and PKS B1039-47 is from 
the compact inner jet regions on parsec scales. This leads us 
to suggest that RM time-variability in extragalactic point 
sources may be due to the evolving polarized jet struc- 
ture on parsec scales which illuminates different parts of 
an inhomogeneous magneto-ionic medium in the immediate 
vicinity of the jet. Follow-up observations of these particular 
sources with parsec-scale spatial resolution using the Aus- 
tralian Long Baseline Array (LBA) will enable us to test 
our predictions. Hence, with multi-epoch polarization ob- 
servations using wide-bandwidth facilities like the ATCA it 
may be possible to map out the parsec-scale evolution of 
polarized components as well as the Faraday rotating envi- 
ronment in AGN jets. 

In the near future, combining data from the pristine RFI 
environment of ASKAP from 0.7 to 1.8 GHz with data from 
the ATCA at higher frequencies can provide an exquisite 
probe of the polarization properties of a much larger sample 
of AGN. 
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Figure 3. Plots of the RM spectra for all four sources. Dashed line: RM synthesis spectrum. Solid line: RMCLEAN spectrum with the 
locations of the clean-components also shown as vertical lines. The text in the top-left corner of each panel lists the parameters extracted 
from the peak in the RMCLEAN spectrum as well as the value of the reference wavelength, Aq. The asterisks connected to the Faraday 
depth axis by the solid line denote the locations of the best-fit model RM. 




0.00 0.01 0.02 0.03 0.04 0.05 0.06 0.07 -0.06 -0.04 -0.02 0.00 0.02 0.04 0.06 

A- [m^] 1 



Figure 4. Polarization data for PKS B1903-802, and the corresponding best-fit single RM-component model. (Eqn. |8)l. Top left: q (open 
circles) and u (full circles) data vs. A^, fitted with the model q (dot-dashed line) and u (dashed line). Top right: p vs. A^ data over-plotted 
by the model (solid line). Bottom left: '3/ vs. A^ data over-plotted by the model (solid line). Bottom right: u vs. q data over-plotted by 
the model (solid line). 
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Figure 5. As for Figure |4] for PKS B1903-802, but now modelled by a single RM-component model with depolarization from external 
Faraday dispersion (Fan. [TTjl . 
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Figure 6. As for Figure[4]for PKS B1903-802, but fit witli a two RM-component model. 
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Figure 7. Polarization data for PKS B0454-810, and the corresponding best-fit single RM-component model (Eqn.[8ll. Layout as described 
in Figure |4] 
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Figure 8. As for Figure [7] for PKS B0454-810, but fit by a single RM-component model with depolarization from external Faraday 
dispersion (Eqn.[TT}. 
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Figure 9. As for FigureElfor PKS B0454-810, but fit with a two RM-componcnt model. The kink in the q(A^) vs. w(A^) plot is a unique 
signature of this particular model. 
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Figure 10. Polarisation data for PKS B1610-771, and the corresponding best-fit single RM-component model (Eqn. [S]!. Layout as 
described in Figure |4] 
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Figure 11. As for Figure [Tol for PKS B160-771, but fitted by a single RM-component model with depolarization from external Faraday 
dispersion (Eqn.[TT}. 
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Figure 12. As for Figure [TOl for PKS B160-771, but fitted with a two RM-component model. 
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Figure 2. Plots of the rotation measure spread function (RMSF). 
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Figure 16. Simulated three component model without noise for 
PKS B1039-47 but with identical coverage as in Figure (S^d) . 
Stars: input model components taken from the best fit to PKS 
B1039-47 data shown in Figure[T4l Dashed line: Dirty RM synthe- 
sis spectrum from simulated data. Solid line: RMclean spectrum 
from simulated data. Clean-component positions and amplitudes 
are also included on plot. 
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